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Abstract 

Semigroups describing the time evolution of open quantum systems 
in finite dimensional spaces have generators of a special form, known as 
Lindblad generators. The simple generators, characterized by only one 
operator, are analyzed. The complete set of all the stationary states is 
presented in detail, including a formula to calculate a stationary state 
from the generating operator. Also the opposite task can be fulfilled, to 
construct an evolution leading to a prescribed stationary state. 
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1 Introduction 

Complete positivity is an essential property of maps acting on states of open quan- 
tum systems, representing their change in time. It was in the 1970s, starting with 
|K71j . when this has been discovered. In the sequel the generators of semigroups 
formed by maps with this property have been identified for finite dimensional 
systems, and also for infinite systems, provided a strong continuity property is 
assumed in addition to complete positivity, |GKS76l IL76] ; they are now known 
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as "Lindblad generators" . The interest at that time has been mainly in connec- 
tion with Non Equihbrium Thermodynamics, see |D76l IS80t lAFOlt IBP02] and 
references therein. 

New interest has arisen recently in connection with quantum engineering of 
small systems, and it seems necessary to enlarge our knowledge about such semi- 
groups. So we proceed in the analysis of these completely positive semigroups, 
acting on states for systems with finite dimensional Hilbert spaces. We are inter- 
ested in characterizing both the set of invariant states and the paths of changing 
states. The early results on invariance stated in \S76\ \F78\ IS80j concern systems 
under special conditions on the set of generators. We give more general detailed 
insight, not needing these restricting conditions. Pioneering work has been done 
by K. Dietz |D03l \^0^ IDOSj . identifying "superselection sectors", discussing the 
change in time of entropy and presenting a formula which relates a stationary 
state to the operator which appears in a simple generator. In the present paper 
his findings are imbedded into a complete mathematical characterization. In ad- 
dition we shed some light on the approach to equilibrium, on the geometry of the 
paths in the set of states. 

There are three types of physical processes one wants to describe: Decay, 
dissipation and decoherence. Decay of excited states leads - in an idealized model 
- to a complete annihilation of some subspaces of the Hilbert space. Dissipation, 
on the other hand, leads to a spreading of the state over a large space, getting 
more mixed, in case it started as a pure state. It is somehow the counterpart 
to decay and also to conservation laws. In reality and also in our studies both 
of these effects occur simultaneously, and conservation laws may impose some 
structuring. Decoherence in a wide sense is dephasing; in a preferred basis the 
off-diagonal elements of the density matrix are diminishing and finally vanishing. 
Each type of these processes has its mathematical expression in properties of the 
quantum dynamical semigroup's generator. 

In this paper we concentrate on simple generators, the building blocks for 
generators of general completely positive semigroups - a characterization of Com- 
plete Positivity is presented in the Appendix. For two dimensional systems, a.k.a. 
qubits, we present a detailed description. In a companion paper we extend our in- 
terest to the general quantum dynamical semigroups on finite dimensional Hilbert 
spaces. 

Our starting point, the connection with the earlier studies, is the result of 
[GKS761IL76] : 

1 PROPOSITION. Generators of semigroups: Every generator of a semi- 
group of completely positive trace preserving maps T* : p{s) ^ p{s + 1) for t > 
on the set of finite dimensional density matrices p, can be written in the form 

p = V{p) = -i[H,p] + Y,Vy,.{p), where H = H\ (1) 

a 

where the dissipative parts, the T>h 
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2 DEFINITION. Simple generators: 

Vh{p) = hph) - \{h^hp + ph^h) (2) 

The commutator with a Hamiltonian alone generates a unitary time evolution, 
which is well known in principle. We concentrate on the other building blocks, 
the simple generators and the evolutions they generate. 

It is well known that the division of T) into a sum of several simple generators is 
not unique. Different sets of operators /i" can be attributed to a given generator 
v. The attribution of a single /i to a simple generator on the other hand is 
mostly unique, with few exceptions. The method of proof for this statement 
uses a structuring of the Hilbert space and of the operator /i, which reflects some 
characteristic properties of the simple V. These methods are presented in Section 
13. H At the end of the same Section the mentioned uniqueness is stated in the 
Lemma [7] and then proved. 

Our analysis starts with studying the decompositions of and their relations 
to the action of the generator. It turns out that the time evolution of p in 
the convex set S of density matrices can be analyzed by studying the action 
of T> on the boundary of S, which is the set of density matrices with rank less 
than the dimension of the Hilbert space. The structure of the operator h and the 
corresponding structuring of the Hilbert space are here of prime importance. Only 
a density matrix the support of which is an eigenspace of h does not dissipate 
in first order of t into higher dimensional subspaces. Therefore we call these 
faces of S which are formed by the p with support being an eigenspace or a 
generalized eigenspace of h the "lazy faces". An even stronger condition marks 
certain eigenspaces and the corresponding faces as enclosures. No path p{t) leads 
into or out of them, some conservation laws are valid. A closer look on the action 
of T>h on density matrices in special faces lets us determine precisely the complete 
set of stable, "stationary", states. These results are presented in Theorem [TH A 
note on naming: In the literature the characterization of a state as "invariant", 
"stable" or "stationary" is used, each one meaning p = 0. In this paper we 
favor "stationary" . And we are sloppy in not differentiating "state" and "density 
matrix" . 

Further analysis, making a detour using the dual time evolution of the oper- 
ators and the Kadison inequality, leads to the insight that there are no periodic 
evolutions, no circular paths for the p{t). All paths in the space S lead to the 
stationary states. The entropy of the states along a path may increase or decrease 
and may even pass through several local maxima and minima. 

Some examples, most of them with qubits, demonstrate the findings of our 
analysis. 

At various places we have to do calculations concerning a splitting of the 
Hilbert space into two orthogonal subspaces. The elaborated formulas are col- 
lected in Appendix 16.21 
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2 The set of generators 

2.1 Basic properties of the generator's actions 

By inspection of the definition ([2]) it is obvious that the semigroups generated by 
the superoperators Vh preserve the trace and selfadjointness of matrices. 

jTrpit) = TT[V{pm = 0, |at(t) = V{a\t)) = [I?(a(t))]t. (3) 

To see that they also preserve positivity, we look at the parts of the simple T>h 
separately. The first part, mapping a positive p to a non-negative operator hph\ 
causes an increase of p. The second part causes a decrease of the eigenvalues. 
The argument goes as in proving the Feynman-Hellmann theorem and gives: 

3 PROPOSITION. A differential inequality: Ifr{t) is a non-negative eigen- 
value of p{t), its change in time is bounded from below as 

m > -whwMt). (4) 

Proof. Differentiating the eigenvalue equation p{t)il){t) = r{t)ilj{t), - Rellich's 
theorem guarantees that all terms are analytic functions, - and using = 

gives 

r{t) = {^\p\^) = {^\{hph^ -^{h^hp + ph^h))\ilj) 
> -l{^\{h^hp + ph^h)\^) 

= -{ilj\h^h\-^)r{t). (5) 

□ 

This implies that the positive eigenvalues Vj of the density operator obey in 
the course of time the inequality 



r 



(t) > exp {-\\h\\h) ■ rj{0) > 0, (6) 



and positive eigenvalues stay positive. By continuity of the time evolution, zero 
eigenvalues can not get negative, the positivity of density matrices is preserved. 

So the evolutions do not leave S. Moreover, there exist stationary states, i.e. 
states which do not change in the course of time. 

4 PROPOSITION. Existence of stationary states. For each simple gen- 
erator there exists a stationary state poo satisfying ^^{p^) = 0. If h has zero 
as an eigenvector, each state with support in the zero subspace is stationary. If h 
does not have zero as an eigenvalue, one may construct a stationary state by 

Poo = {h^h)-^/i:v[{h^h)-\ (7) 
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Proof. If p = X]j "where h\(j)j) = 0, obviously hp = 0, p/i"^ = and 

therefore 'Dh{p) = hold. For h without zero as an eigenvector, consider the 
polar decomposition h = U\h\, where \h\ = y/Wh and U = h\h\~^. The evolution 
equation for poo = |^|~^/Tr[|/i|~^] is 

Poo = U\h\p^\h\U^ - i(|/i|Voo + Poo\h\') = {U-U^- l)/Tr[|/i|-2] = 0. 

□ 

This relation between h and poo has already been given in |D03j . The questions 
of uniqueness and of attraction to paths of other states are answered in Theorems 
US] and [13 

2.2 Superoperators 

In the sequel it will be useful to endow the linear space of matrices with an 
inner product (o"|p) = Trcr''"p, and consider matrices as elements of this space, 
the Hubert- Schmidt, HS space. In the middle of the set of density matrices S 
is the maximally mixed state u = 1/ dim(7i). The square distance ||p — a;|p is 
related to the "Fermi entropy" Trp(l — p), on which we make some remarks in 
Subsection 14.31 

Consider a generator D as a superoperator acting on the HS space. In this 
respect the semigroup which it generates can be extended to a group. For neg- 
ative t the trace and selfadjointness of matrices are still preserved, but not the 
positivity. 

The adjoint V'^ can be looked upon as generator of the dual time evolution of 
observables in the Heisenberg picture. 

Vi{f) = h^fh-\{h^hf + fh^h). (8) 

Preservation of the trace of density operators corresponds to preservation of the 
unity under the dual time evolution of observables, the semigroup generated by 
pt is "unital" . (That both V and 

are unital and trace preserving occurs, if h 

is a normal operator.) 

Diagonalizability of V is rather exceptional. Nevertheless there exists for 
each D its set of eigenvalues {A}, equal to the set of zeros of the characteristic 
polynomial det(T' — Al), a set of proper eigenmatrices {a}, obeying Va = Act, 
and, in some cases of degeneracy, generalized eigenmatrices, obeying ("D — A)"cr = 
0. As in any set of linear differential equations, any p(t) can be expanded into 
a sum of such a ■ e^^ ■ Polynomial (t). Some general properties of eigenvalues and 
eigenmatrices are presented in Subsection 13. 4[ Examples and special cases will 
be presented for twodimensional systems in Subsection 14.21 

The generators V, defined in Proposition [T], with H = form a convex cone. 
The simple generators make the extremal rays. Nevertheless closure adds the 
generators of unitary groups: 
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5 PROPOSITION. Generators of unitary time evolution can arise as limit of 
simple Lindbladian generators. 

Proof. Define h{X) := 2X~^1 — iXH where = H. Consider the limit A — 0: 
Vhix)ip) = -t[H,p] + X'Vnip) ^ -t[H,p] (9) 

□ 

Note that the norm of h{X) diverges, but 'D^x) converges in norm. 

So it is to be expected that a thorough inspection of the extremal rays, con- 
taining the building blocks, the simple generators, will give a good basis for the 
analysis of all GKS-Lindblad generators. 

3 Simple generators 

3.1 Splitting, Decomposition, Structure 

We consider the generators V = V^, defined in Definition [21 To analyze their 
properties it is convenient to find structures on two levels. The structure in the 
large concerns a kind of enclosure: subspaces of Ti. with the property that the 
system can neither leave nor enter any of them in the course of time. (These 
subspaces are related to the superselection sectors in |D04j .) P is a projector 
onto such a subspace if and only if it obeys 

[P,h]=0. (10) 

Each operator h appearing in a generator Vh can be decomposed as 

h = ^h, (11) 
j 

with a decomposition of the Hilbert space Ti. = @j Tij into mutually orthogonal 
subspaces Tij = PjH, where each hj acts on Hj and is not further decomposable. 

Aim and purpose of this decomposition in the large is to reduce the analysis 
of the superoperator's action on any p to the analysis to the split action on the 
blocks PjpPj and PjpPi. It is: 

I), {PjpPj) = h,ph] - \[h]h,p - ph]hj\ = (P^pP,) , (12) 

for the diagonal blocks in the decomposition of p into block matrices PjpPe- The 
off-diagonal blocks give the phase relations between the enclosure-like subspaces. 

Vh (PjpPe) = hjphl - i [h]h^p - ph\he] . (13) 
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In case of an appearance of equivalent parts hj there is the possibihty of 
unitary reshuffling the corresponding subspaces: Collect the set of equivalent 
parts as /ii . . . /iat, i.e. 



Represent their sum hi® ■ ■ - (BfiN as a tensor product 1 at ® /ii acting in C ® ?i i . 
The action of a unitary operator V © 1 in this tensorial representation gives a 
unitary reshuffling of eigenspaces Hj. 

6 PROPOSITION. Splitting of operators and of spaces: The splitting 
h = hj with a set of orthogonal projectors Pj which commute with h is unique, 
except for a unitary reshuffling of eigenspaces Tij = Pj'H on which equivalent parts 
hj are acting. 

The action of the generator splits as 



Proof. By induction, one may find a maximal set of mutually commuting orthog- 
onal projectors P,, all commuting with h: If there exists one P commuting with 
/i, one may start with the set {P, 1 — P}. Then one searches for a further splitting 
P = P2 + P3, the Pj commuting with h. If such a further splitting exists, one 
replaces P by P2 and P3. The same for 1 — P, and then further splittings of the 
Pj, and so on. At the end of this procedure, one has the splitting Ti = 0^- PjH, 
and h = PjhPj, which can not further be refined. 

Now collect piecewise all the mutually equivalent parts /ij^, as l„(j)®/ij^ 

acting in C"'--'^ AH the possible families {Pj} of projectors can be observed 

in this way: The operator h commutes with each projector onto V (S> where 
V is any subspace of C"'^^\ and h commutes with no other projector. 

The split action of the parts on diagonal and off-diagonal blocks of p, stated 
in the equations f|T2l) and f|T3l) . is obvious. □ 

This decomposition is an important step in the analysis of the time evolution. 
An algebraists point of view is that the projectors commuting with h create 
a subalgebra of conserved observables, T>'^{Pj) = 0. The law of conservation 
implies that the probability for the system to remain in an enclosing subspace 
PjTi. is constant in time. The evolution inside the subspace is determined by the 
action of . Studying the total action of Dh consists of determining the action 
of single Vh. with indecomposable hj, and studying the evolution of the phase 
relations between different subspaces Tij. The indecomposable parts acting on 
diagonal blocks of p are studied in Sections [3.21 and [33| The evolution of the phase 
relations, represented in the off-diagonal blocks PjpPe, is discussed in Sections 



The structure on the second level is related to representions of an indecompos- 
able h - which may enter in the large as an indecomposable part hj - in a Schur 



V{j, i} C {1...N} 3 unitary Uje : hj = UjehU] 




(14) 



EJlandim 
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triangulated form, as an upper triangular matrix, see f.e. |L69] : On the diago- 
nal are the eigenvalues in arbitrary sequence, below the diagonal there are only 
zeros. An indecomposable h of dimension greater than 1 is not diagonalizable, 
one has to be aware of the generalized concepts: The eigenspaces and general- 
ized eigenspaces are subspaces PTC with projectors P, with the characterizing 
property 

hP = PhP. (15) 

These P are orthogonal projectors, P = P^; each one-dimensional P projects onto 
a proper eigenvector. These eigenvectors are in general not mutually orthogonal. 
The eigenspaces and generalized eigenspaces are spanned by the proper and the 
generalized eigenvectors, ip\j,o and ip\,j,m, of h. 

{h - = 0, {h- X)lp\,j,m = ^A,i,m-l 

Note that a generalized eigenspace containing ipxj^m has to contain also all the 
V'Aj.^i with fi < m. 

To be sure about the relevance of the properties of h for the action of Vh it 
is helpful to study the uniqueness of attributing h to a given simple generator. 
Already here the method of Schur triangulation appears as an important tool. 

7 LEMMA. If it is possible to write a generator D as a simple T>h, the attribution 
of h is unique up to a phase factor, unless h is a "normal", i.e. diagonalizable, 
operator. 

Proof. The strategy of the proof is to consider the equations p = T>h{p) for the 
given h, with a set of properly chosen matrices p as determining the operator k 
which shall give the same p, by p = Vk{p). The condition of non-diagonalizability 
of the given h can hold only for n = dim{T-C) > 2. The proper choice of p requires 
a proper choice of the basis vectors, giving a triangular matrix representation 
of h and enabling a proof by induction on the rank of p. If the given h has a 
generalized eigenvector ipx,!^ choose the basis for its decomposition and its Schur 
triangulation in such a way that ipx^ is the first basis vector and ipx^i is the second 
one. If h has only proper eigenvectors, then, because of its non-diagonalizability, 
there must be a non-orthogonal pair of them, ipxfl and ^/'^.o- Choose this pair 
to form the first two basis vectors: ipxfi and the normalized tp^^Q — 'Jipx,o, where 
7 = {'^x,o\'^fj,,o)- Then we proceed by induction on the dimension of the space and 
use the chosen basis. We refer to the formulas presented in the Appendix 16.21 

(r Q \ 
+ with (7 = 

s J 

and s = determines c = for the matrix k = ( ^ j ] • The action on p with 



^ c d 

r = q = determines first which is not zero, due to the proper choice of the 
basis vectors. Then via q = a^br it determines a^b. Finally, acting on p with 
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r = q = 0, hj q = {hd) — ^a^h)s it determines hS . So the matrix k is determined 
up to a phase as k = h. 

For higher dimensions use the block matrix form, as in the Appendix 16.21 with 
rank(74) = rank(i?) = n — 1, where S = s and D = d are numbers, B = \b) and 
Q = \q) are vectors. The action on a p with \q) = and s = is R = Va{R) 
and - there is the induction hypothesis - this determines A up to a phase. Note 
that R 0, so A ^ al. The action on this p is also s = 0, determining C = 
and \q) = RA^\b). Choosing a phase for A, the action on p with i? = 0, |g) = 
determines \b) via R = \b)s{b\ and via \q) from before; and it determines d via 
\q) = s{d^b) — ^A^b)) in case \b) ^ 0. If |6) = 0, consider the action on a p with 
|g) 7^ 0: \q) = d^ A\q) — }^J^ A + |c?p]|g). This is a quadratic equation for the 
number (i, but there are at least two different equations for the components of |g) 
to determine d. So the proof of uniqueness is completed for nondiagonalizable h. 

□ 

\i h\s a diagonal operator, all the matrix elements of p evolve in time inde- 
pendently of each other, 

= - Kl^iP + l^jHlPij, 

the diagonal of p is invariant in time. Only diagonal k in can give the same 
evolution. But there are not always enough equations to determine the matrix 
elements on the diagonal of k. There remains some freedom in choosing k with 
\ki — kj\ = \hi — hj\ and Im{kik'^-) = Im{hihj) in special cases, e.g. if there are 
not more than two different values of hj, or if all hj are real. But the qualitative 
structure of possible k is fixed: It is determined that it is a diagonal operator in 
the same basis, with the same degeneracies. 

3.2 At the boundary of the set of states 

It turns out that the time evolution of p in the set S can be analysed by studying 
the action of V on the boundary, the set of density matrices with rank less than 
the dimension of the Hilbert space. Referring to the inequality ([6]), we note that 
it implies one of the central observations: 

8 LEMMA. No purifying in finite time: The rank of p{t) cannot decrease. 

An observation of almost sure increase of the rank serves as a companion to 
this exclusion of a decrease. 

9 LEMMA. Consider p{t) with rank(p(0)) < dim(7-^). At least one eigenvalue 
which is zero for p(0) becomes strictly positive in first order in t, if and only if 
the range of p{fS) is not an eigenspace or a generalized eigenspace of h. 
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Proof. Consider the analytic expansion of p{t) and its eigenvalues r{t), using the 
eigenvectors related to them, and consider the terms linear in t. 

Write h in Schur triangulated form, such that the upper left part R (we use here 
the notation of the Appendix 16. 2p acts as a matrix of full rank in the range of 
p(0). Use the block matrix equation (HHl) . with R the restriction of p to its range, 
an eigenvector to the eigenvalue zero, so it is not in the range of R and 

f(0) = (01^(0)10) = (0|Ci?Ct|0). (16) 

Observe that, according to (ITSl) . the range of p(0) is a (generalized) eigenspace, 
iff C = 0. Since R is of full rank there is at least one giving a strict increase 
r > 0, iff C ^ 0. □ 

Note the correspondence between subspaces of the Hilbert space Ti. and faces 
of the set of states S: Each face consists of density matrices whose range is some 
subspace of H. Since not all subspaces of H can be eigenspaces of a non-constant 
operator, the corresponding faces form a zero-set in the boundary of the set of 
states. Here comes the structure of the operator h into play. 

10 DEFINITION. We define each face which corresponds to an eigenspace or 
a generalized eigenspace of h as a lazy face. In a formal notation: If P is the 
projector onto the range of p, and if P fulfills the property fTS\) i. e. hP = PhP, 
then, and only then, is p in a lazy face. 

The reason for this definition is Lemma [9], and its strengthening in the follow- 
ing Lemma [m Eigenvalues r{t) which are zero for t = do not increase in first, 
not even in second order in t if p(0) is in a lazy face. They may increase in third 
or a higher order, or not at all. Examples will be presented in Subsection 14. 2[ 

11 LEMMA. Consider p{t) with p(0) in a lazy face. Each eigenvalue r{t) be- 
ginning as r(0) = does not increase in first and not in second order in t. 

Proof. The eigenvalue equation p{t)il){t) = r{t)%l){t) has to hold in each order of t. 
Consider the splitting = x©0, with x(0) = 0, the block matrix form presented 
in Appendix 16.21 with C = 0, Q(0) = 0, S'(O) = 0, and the Taylor expansions. 
For the blocks of p we need 

R{t)=R + 0{t), Q{t) = -\RA^Bt + 0{f), S{t) = ARA^ Bt^ + 0(f). 

(17) 

The first order of the eigenvalue equation, p(O)'0(O) + p(O)-0(O) = r(O)'0(O), gives 
two equations, one for the vector component orthogonal to the range of p(0). 



-\B^ARx{Q) = r(O)0(O) r(0) = 0, 
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and one for the component in the range of p(0) 

-lRA^B(l){0) + Rx{0) = x(0) = \A^B(j){0). 

The last imphcation uses the invertibihty of R, due to its full rank. From the 
second-order equation we need only the part in the subspace orthogonal to the 
range of p(0): ^p(0)?/'(0) + p(0)V^(0) + |p(0)V^(0) = |r(0)V(0) gives 

^B^ARA^B(f){0) - iB^ARxiO) + = |r(O)0(O) f(0) = 0, 

since the terms on the left hand side cancel when the identity for x(0) is used. □ 

There is a way to geometrically characterize the action of V on the states in 
a face. Define the mid point Uf oi the face as 

j 

where the (pj form a basis of the n-dimensional subspace corresponding to the 
face. Consider the projection of the path p{t) in the Euclidean space of hermitian 
matrices onto the line which connects uj with the maximally mixed state u in the 
center of S. This projection gives a special coordinate for the HS vector pointing 
from Uf to p{t): 

TT[ip{t)-ujf)iuj-ujf)] = ^TT[S{t)]. (18) 

Only for lazy faces does this coordinate not grow linearly in t. It grows quadrat- 
ically in t, see ( !T7|) . so this geometrical characterization is not as strong as the 
statements on the eigenvalues made in Lemma [TTJ 

In the next Section we explain in detail how an appearance of zero as an 
eigenvalue for h has consequences in the investigations on the set of invariant 
states. They can be elements of the corresponding lazy face. Concerning the 
dimension of such a face we state as a remark the following 

12 LEMMA. An indecomposable operator h of rank n can not have more than 
n/2 proper eigenvectors to a degenerate eigenvalue. 

Proof. Write h in triangulated form and use the m proper eigenvectors to the 
degenerate eigenvalue A as the first m basisvectors. 

d)' 

where A = Xlm- The rectangular matrix B can not have more rows than columns, 
i.e. m < n — m; otherwise the rows could not be linearly independent and there 
would be a unitary transformation of the m-dimensional proper eigenspace such 
that m — (n — m) rows of the transformed B would be zero, contradicting the 
undecomposability of h. □ 
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3.3 Stationary states 

The upshot of the investigations on the action of at the border are the following 
two complementary propositions: 

13 PROPOSITION. Quitting the boundary. has no stationary state at 
the boundary if h is indecomposable and does not have zero as an eigenvalue. 

Proof. Consider p at the boundary, write it as a block matrix (notation as in the 
Appendix 16. 2p with R a matrix of full rank, i.e. R has no zero-eigenvalues, Q = 
and S = 0. If C 7^ 0, use equation (H8|) . which implies 5* 7^ 0. Otherwise, apply 
equation fl50l) : Q = RA'^B 7^ 0, since neither R nor A'^ have zero eigenvalues and 
-B 7^ because of the non-splitting property of h. □ 

But otherwise, one of the lazy faces may be an attractive face: 

14 PROPOSITION. Eigenvalues zero are attractive. // h is indecompos- 
able and has zero as an eigenvalue, then the proper eigenspace to this eigenvalue 
zero corresponds to a face with stationary states. This face is an attractor, and 
no other state is stationary. 

Proof. If p is in this face, obviously hp = 0, p/i^ = and therefore T>h{p) = 
hold. On the other hand, if p is not in this face, one analyzes P'^pP'^, where 
P is the projector onto the proper eigenspace belonging to the eigenvalue zero 
of h. Use the triangulated block matrix form presented in Appendix 16.21 with 
A = acting on the zero-subspace, with C = 0, and S = P-^pP^. Equation (|5Til 
reduces to S = Vd{S)-\[B^ BS + SB^ B], and |Tr[5] = -TilBSB^]. Looking at 
5*00 = lim-r^oo -f Jq S{t)dt, which must be stationary in time, we find that both 
terms have to vanish separately in 6*00 = ^£,{800) — \ [B^ BSco + SooB"^ B] = 0. 

First, assume h does not have improper eigenvectors to the eigenvalue zero. 
This means that all the zeros in the diagonal of the triangulated h appear in A 
and D = P^hP^ is without zero-eigenvalues. Since h is indecomposable, i? 7^ 0. 
If p has S* 7^ 0, then, as is stated in proposition [T31 V^^S) 7^ unless S is 
a matrix of full rank without zero eig envalues; but then B^BS + SB^B ^ 0, 
acting as emptying out the subspace which is orthogonal to the zero-subspace, 
and leading to S" = 0. 

Now it remains to consider h with improper eigenvectors to the eigenvalue 
zero. These zeros appear in the diagonal of D = P^hP^, at least one of them as 
a proper eigenvalue of D, which acts in P-^Ti, a space of lower dimension than 
that where h acts. We make an induction on the dimension of the Hilbert space. 
The induction hypothesis is, that the sub-subspace /C spanned by the proper 
eigenvectors for the zeros of D makes an attractive face for S{t), which is defined 
by -^S = one part of the time evolution of S{t). The stationary 6*00 

vanishes on 7i — /C. But the sub-subspace K. is an improper eigenspace for h, for 
each \^) e /C one has B\ij) ^ 0. So |Tr[5(t)] = -Tr[BS{t)B^ ^ for each S 
which does not vanish on /C, and P-^Ti. is emptied out in the course of time. □ 
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15 THEOREM. All the stationary states for simple generators. 

a) Dissipation. If h is indecomposable and has no zero-eigenvalues, there is a 

unique stationary state in the interior. It is 

Poo = {h^h)-'/TT[{h^hr'] (19) 

b) Decay. If an indecomposable h does have zero- eigenvalues, the stationary 

states form an attractive face, corresponding to the space spanned by the 
zero- eigenvalues. 

c) Elementary dephasing. If, in the other extreme, h is diagonalizable, the 

set of invariant states is the set of density matrices commuting with h. 

d) Stationary splitting. If h is decomposable, h = ^jhj, each hj indecom- 

posable, the direct sums of stationary states pj of its parts form the set 
of all the stationary states which do not have phase relations between the 
independent subspaces Hj. 

e) Dephasing of parts. If h is decomposable, there is a set of invariant 

phase relations pj^i between the subspaces which are domains of definition 
of hj and hi, if and only if either both hj and hi have zero-eigenvalues, 
or hj = hi. In the first case the invariant off-diagonal block matrices are 
Pj/ ~ Ylioi i3^oL,i3\vj,o){wi^i3\ where the \vj^a) cind jw^,/?) are the proper zero- 
eigenvectors of hj and hi. In the second case any invariant phase relation 
Pj^i is equivalent to pj multiplied with some complex number, where pj = pi 
is the stationary state ofT>hy 

Proof, a) The stationarity of poo has been stated and proved in Proposition [Hand 
in its proof. Suppose there were two different stationary states, po and poo- By 
linearity of the evolution equation, the whole line po + A(poo — Po) would consist 
of invariant matrices. But this line would intersect the boundary where there are 
no stationary states, see Proposition [T^ 

b) See Proposition [m 

c) That [/i,p] = implies 'P/i(p) = is obvious. The off-diagonal elements 
obey pij = {h*hj — \{\hi\^ + |/;,jp)), which is not zero, if hi ^ hj. For details on 
the decrease of phase relations pij in this case see Subsection 14. 1[ 

d) Thats the combination of the previous results, using Proposition El 

e) Any density matrix p for the whole Hilbert space Ti can be considered as 
a big block matrix with entries pj £, corresponding to the splitting oiTi = 0^. Tij 
defined in Proposition [61 Each pj i has a time evolution on its own, without 
any mixing with the other blocks. The stationary diagonal elements pj^ have 
been identified in a) to d). Knowledge about them, together with the general 
preservation of positivity, helps to find possible invariant off-diagonal blocks pj^i. 
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So, for any pair hj, hi of indecomposable parts of h, we look at a state p restricted 
to the subspace Tij © TY^. 

First assume that hj has zero as an eigenvalue. Since positivity is preserved 
as t — > cxD, and since stationary pjj reside on the subspace spanned by zero- 
eigenvectors, an invariant pj^i must also obey hjPj^ = 0. Its time evolution is 
Pj^£ = —-^pj^ihlhi, and this can vanish only if Pj^ehe = 0. Invariant phase relations 
exist therefore only if also he does have zero eigenvectors. The stationary pj^i can 
then be expanded in ket-bra products of zero-eigenvectors. 

Now consider the case that both hj and hi are indecomposable and do not 
have zero as an eigenvalue. Assume that an invariant pj^i ^ exists. Keeping 
the diagonal blocks pjj and pi^i fixed, given by the unique stationary states 
determined by hj and hg, multiplied with any non-zero weights, the off-diagonal 
elements can be chosen with some factor A, as Xpj^g and pej = X*Pj/, such that 
the whole p is at the boundary of the set of states. Now we refer to the formulas 
stated in Appendix 16.21 formula (HHi) and the following lines. As is shown there, 
the density matrix, now represented by R acting on a subspace, is invariant only 
if C = and also RA^B = 0. Now R has no zero-eigenvalues, and A has a subset 
of the whole set of diagonal elements of hj © hi as eigenvalues, also without zeros. 

So -B = must hold. But then hj(Bhe = ^ ^ commutes with the projector 



1 




Since neither hj nor hi commute on their domains of definition 



PjH and PiH with anything else but the constant operators, the only possibility 
for such a nontrivial P is the unitary equivalence of hj with hf. hj = Vh(V\ 
V -V^ = 1. 



P = plj + Vpil-p)V + ^/p{l-p)V^ + (1 - p)U 

This unitary equivalence implies then the equivalence of the invariant blocks with 
some factors, equivalence of Cj ■ pjj with Cj^i ■ pjj and q • pi^i. □ 

We may reformulate the results, so that equivalent parts are collected as 
factors in a tensor product: If one knows the stationary states pj^a for an inde- 
composable hj, one knows all possible states for l®hj; they are Yla ® Pj,a, for 
any positive matrix cTq, which has Tr[crQ] = 0. Stationary phase relations between 
the parts may exist. 

In the case of an indecomposable h with no zero-eigenvalues it can be shown, 
as a companion to PropositionfT^ that the unique stationary state is an attractor. 
See Subsection 13. 4[ 

Parts a) and b) of Theorem [15] and the formula (fT9l) allow for an inversion, 
constructing an operator h in such a way that the evolution generated by Vh has 
to lead to a given state p. 

16 THEOREM. Construction of a simple evolution leading to a given 
state. For p either a pure state or a state with a density matrix of full rank, 
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but non constant, there exist dynamical semigroups with a simple GKS-Lindhlad 
generator which have the given p as the unique stationary state. 

Proof. In case p is a pure state, take an indecomposable h with just one zero 
eigenvalue, so that hp = 0. This may be, for example, h = J2j I0j)(0j+i|7 with 
|0o)(0o| = P- 

In case p is a matrix of full rank, construct h = U\h\ with \h\ = p~^^'^, 
as it has to be according to formula (IT^ . This h is decomposable iff there 
exists a nontrivial projector P commuting with h. This commutation implies 
also the commutation of P with h'^ , with = h'^h and with \h\. Since this 
\h\ has no zero eigenvalue and is invertible, P has also to commute with U. 
As commuting operators, P and \h\ have common eigenvectors (pj, such that 

1 /2 

\h\ = E,r-^'10,)(0,|, and P = Now one has to consider a 

unitary operator U which mixes all the eigenspaces of \h\ belonging to different 
eigenvalues. This may be, for example, U = \(f)j){(f)j+i\ + |0Ar)(0o|, where N is 
the dimension of the Hilbert space. Such a U does not commute with any P of 
the sort which commutes with h, and h is indecomposable. □ 

For exceptional p, not pure but not invertible, and also for the maximally 
mixed uj, it is possible to find a semigroup with this p as a stationary state, 
however it will not be the unique one. 

3.4 The geometry of paths 

The geometric characterization of the paths p{t) is equivalent to the algebraic 
characterization of the superoperator V, determining its eigenvalues, its proper 
eigenspaces and also the generalized eigenspaces. Logical arguments can work in 
both directions, from geometry to algebra or vice versa. 

On the geometrical side we know that positivity of p is preserved, that the 
compact set S of states is mapped into itself. The algebraic consequence for V is 
that there is no eigenvalue with strictly positive real part. The detailed argument 
for this fact can be seen in the following discussion which is needed for deeper 
analysis. 

The reflection cr ^ o""^ is compatible with the action of T>, i.e. I'(o"^) = [^(o")]^. 
The consequences are: If an eigenvalue A of P is a real number, the accompanying 
eigenmatrix and generalized eigenmatrices - if there exist any - can be chosen as 
selfadjoint. Each pair of a complex eigenvalue A together with the eigenmatrix a 
has a mirror companion in the pair A* and a^. For the set of self adjoint matrices 
this implies that for each complex number z the matrix Tz = za + z*a'^ lies on a 
path Tz{t) which is a spiral. The imaginary part of A determines the time T of 
revolution, the real part gives the change of the HS-norm. 



T := 27r//m(A), r,(T) = e^'^^'^^^^^r. 



(20) 
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Now consider any p in the interior of S, and some number e, such that 



p± = p± er, 



z 



(21) 



are both inside of S. Their distance changes exponentially in time, 



p+(t)-p_(t)|| = ||p+(0)-p_(0)||-e' 



,27r_Re(A)f 



But both of them remain in S, their distance can not increase, Re{X) < 0. 

It remains to discuss the cases of Re{X) = 0. We know for sure that stationary 
states exist. They are eigenmatrices to the eigenvalue A = 0. We state the absence 
of other cases as 

17 THEOREM. No circular paths. There are no circular paths in the HS- 
space of matrices; there are no eigenvalues X of V on the imaginary axis except 
A = 0. 

Proof. We decompose the Hilbert space Ti = @j the operator h = 

PjhPj into indecomposable parts, as stated in Proposition [6l For matrices a 
we use the block matrix decomposition a = iPjC^Pe and observe the mutual 
independence of evolutions for the blocks, V{PjaPi) = PjV{a)Pi. So the search 
for eigenmatrices can be done by considering single block matrices. 

First we consider hj having zero as an eigenvalue. We know by Proposition 
[13] that all the states PjpPj move to the attractive face, where nothing changes 
in time any more. Now consider the path p+(t) starting with po + er^, as above 
in equation (121]) . po(^) moves to a poo on the attractive face, as t ^ oo. Also 
p+{t) has to stay inside of S, moving to the attractive face, where there can be 
no circular path. So also p+(t) has to approach poo, so Re{X) < 0, strictly. The 
path may have the structure of a tornado. 

For a matrix cr = PjpPi, giving the phase relation between Tij and Tie, we 
know that it has to move to the subspace of matrices obeying hja = 0, because 
of preserving positivity for states on the domain Tij © Tie. For such a a the time 
evolution is dictated as d = —^ahlhi. This gives either exponential decrease in 
time, or invariance, or a combination of both, but no circular path, no nonzero 
eigenvalue of V on the imaginary axis. 

Now it remains to consider those a = PjcrPj and a = PjcrPe where neither hj 
nor he do have zero as an eigenvalue. Here we turn to study the adjoint It 
determines the time evolution in the Heisenberg picture. Its set of eigenvalues is 
the complex conjugated set of eigenvalues of V, so these sets are actually identical. 

We use the Kadison inequality, [K52] (See the Appendix 16. ip . 



(22) 



which holds for the mapping F i-^ ^t{F) of time evolution in the Heisenberg 
picture. 
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Consider an eigenmatrix F = PiF Pj, with V^{F) = XF. The adjoint is an 
eigenmatrix with eigenvalue A* Assume X = ir with r G R. Then = 
{e-'''^F^){e'''^F) = F^F. The matrix F^F maps Hj into Hj and is positive. Take 
the projector onto an eigenvector with the largest eigenvalue, which is as 
a density matrix p, and switch for the moment between the Schrodinger and the 
Heisenberg pictures: 

Tr[p{t)F^F] = TT[p<^t{F^F)] > Tr[pF^F] = (23) 

The mean value over the time > has to give the stationary state p on Tij, which 
is unique and has no zero eigenvectors in Tij , see the Theorem [151 

Tr[pF^F] = hm ^ rTr[p(t) F^F] > \\F^F\\. (24) 

On the other hand, obviously TrfpF^F] < \\F^F\\, so Tr[pFtF] = In 
a basis where p is diagonal, the matrix representing F'^F has therefore all diag- 
onal elements equal to its norm, which implies vanishing of all the off-diagonal 
elements, so F^F = \\F^F\\lj. 

li £ ^ j one finds also F F^ = \\F F^li by doing analogous investigations. 
A first consequence is, that the existence of such an eigenmatrix is possible only 
if Hj and He have equal dimension. With F F^ F = F \\F^ F\\ = F we 

conclude II Ft F|| = ||FFt||, so 1/ := F/v/||FtF|| is a unitary operator if Hj = Hi] 
and it is an isometry between Hj and Hi if these are different subspaces. We 
turn back to the Heisenberg picture and multiply the evolution equation ([8]) for 
V from the left by V"^: 

ytpt(y) = V^hlVhj - \V^h\hiV - \V^Vh]hj = V^XV = irlj . (25) 

We insert V between hi and hi, define hi = V'^hiV, take the trace and read 
it as an equation for inner products in the HS-space: 

{hi\hj) = ^{{hi\hi) + {hj\hj)) +irn, (26) 

where n = dim(7ij). For j = £ we can conclude immediately that X = ir = 0. 
For j ^ i this follows from the Cauchy Schwarz inequality and the inequality 
between the geometric and the arithmetic mean 

|(/i,|/i,)|<||/i,||||/i,||<i(||/i,f +||/i,f). (27) 

□ 



4 Special systems and examples 
4.1 Dephasing 

The stable phase relations between indecomposable parts have been identified 
in Theorem (fT5|) . In some special cases we can state details on the integrated 
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evolution: One may restrict the study to indecomposable parts. We use the same 
notation as in Appendix 16.21 applied to the operators in the subspace Hj © He, 
writing A = hj and D = hf for the parts of h. The phase relations are expressed 
in Q = pj^£. They evolve, according to equation (l50l) as 

Q = AQD^-\{A^AQ + QD^D), (28) 

a) Looking at a diagonalizable h, with A = a and D = d just numbers, the 
phase relations, the off-diagonal elements of p, evolve according to equation (|28ll 
as 

^q = -^(\a-d\'^ + {a*d-ad*))q, (29) 

giving exponential decrease, \q(t)\ = exp(— |a — d\'^)q{0), together with a phase 
rotation. 

b) Consider D = hp having zero-eigenvectors and Q-D^ = 0. Then 

Q = ~\A^AQ. 

Expand Q = CaQa with A^AQa = daQa- Exponential decrease is the conse- 
quence: 

a 

c) Consider as one- dimensional, with hf = d, and the other part hj = A 
as indecomposable. Then equation (125]) implies monotone decrease of Tt{Q''Q): 

jJ^Q^Q) = -Tr(gt[(At - d^){A - d)]Q) < 0. (31) 

This inequality is strict, unless d is an eigenvalue of A and the columns of Q are 
corresponding eigenvectors. In case of {A — d)Q = this first order derivative of 
HQ IP vanishes. Nevertheless one gets 

Q = (|d|2 _ - lldHQ = ld{d* - A^)Q ^ 0, (32) 

since A is an indecomposable hj, so A and A^ have no common eigenvalues. There 
is no stationary phase relation. The second order derivative of Q^Q in t is also 
zero, but the third order derivative is strictly negative. 

To make the calculation transparent, we consider Q as an element \Q) of 
the Hilbert-Schmidt-space, with the inner product {P\Q) = Tr(P^(5). Time- 
derivative is a super-operator V, mapping \Q) to \Q). In the formulas we use the 
same letter A for the special super-operator which we define as A\Q) := \AQ), 
as for the operator in H. 



\Q) = V\Q) = [-1{A^ - d*){A -d) + '^{Ad* - AU)]\Q). 



(33) 
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Its adjoint operator is 

(Ql = {Q\V^ = (g|[-|(At - d*){A -d)- \{ad* - A^d)]. (34) 
In this formalism we get ^TrQ+Q = {Q\{V^ + V)\Q) and 

— TrgtQ = 5^ M {Q\V^'^V--\Q) (35) 

m=0 ^ ^ 

The assumption = dQ gives (r't+i))|Q) = o and its adjoint, {Q\{V^+V) = 0. 
Now we write the second and first derivatives, using fl35l) . as 

^Trgtg = {Ql ((I?t + + pt(pt + p)) |g) = g, 

and, using |g)(g|/||g||2<i, 

— Trgtg = {Q\ (^(Pt + p)i)2 ^ pt(pt + + (i)t + |g) 

= {Q\V\V^ + V)V\Q) = -{Q\V\A^ - d*){A - d)V\Q) 

< -{Q\V\A^ - d*)\Q){Q\{A - d)v\Q)/\\Qr = -^did* - At)g||7||g|p < 0, 

with a strict inequality, as in equation (l32i) . 

We warn of thinking to generalize these examples. There is in general no 
monotone decrease of phase relations, not in the form of a decrease of the HS- 
norm of Q. Some remarks are stated in Subsection I4.3[ 



4.2 Twodimensional systems 

We present the stationary states, the investigations on the eigenvalues and eigen- 
matrices of T>, and make remarks on the paths. 

If h is constant, then Vh is zero. This case is nevertheless worth mentioning. 
In the way of decomposing a larger system it may appear as a restriction to a 
subspace. The constant h is then active in the evolution of the phase relations 
with other subspaces. 

For non constant h consider the triangulated representation 



• If b=0 there are four eigenmatrices of V^, 

\0 J ' \0 1 J ' \0 J ' \l J 
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the first two belonging to the eigenvalue 0, the other two to A = ^|a — + 
^{a*d — ad*) and to A*. Each diagonal p is constant in time, dephasing occurs as 
stated in Subsection 14.11 part a. 

• For any h with b 7^ there is a unique normed stationary state p, with 

lap -a*b 



P 



r q 



q* S J' ' |a|2+ |5|2 + |^|2 ' - |a|2+ |5|2 + |^|2 

(37) 

This is easily seen by considering the equations ( H9|50f51l) in the simpler form with 
numbers instead of matrix-blocks; and it is, of course, consistent with equation 

Any kind of p can appear as a stationary state for some Vh, except the 
maximally mixed u = 1/2. The stationary po is a pure state if a = or = 0. 

r q 



Concerning all the eigenvalues A of Vh and its eigenmatrices a 

' p s 

we note that the preservation of the trace takes the form s+r = A(r + s) = 0. The 
case A = Aq = involves the unique stationary state. The other eigenmatrices 
belong to Aj 7^ and have r + s = 0. By using the equations ( l50|52|l5ip . together 
with this condition Tro" = 0, one gets the eigenvalue equation & = Vh^a) = Act 
in the form 



dt 



ad*-\\h\\l/2 

a*d-\\h\\l/2 













p 








{ s ) 







(38) 

The remaining eigenvalues of P are the zeros of the characteristic polynomial for 
the 3x3 matrix appearing in this equation. 

• For a=0 there are only real eigenvalues. They are Ai = — A2 = A3 = 

— II /i II 2/2, the eigenmatrices are 

/ \b\'-\d\' 2bd* \ 1 \ \ 



2b*d \d\'-\b\' 

For a = and |6| = \d\ ^ all three zeros of the characteristic polynomial 
coincide, and the eigenmatrix ai is no longer linearly independent of a2 and a-^. 
But there is now a generalized eigenmatrix 



0-1 



1 
-1 



satisfying the equation V^cxi) + \b\^o-i = |6p(o-2 + (T3). 

• For a^O,b^O, d^O (note that c? = is equivalent to a = 0) the three 
eigenvalues of V which are not zero have to be calculated by finding the zeros of 
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the characteristic polynomial for the matrix appearing in (1381) . a polynomial of 
third order. A simple case occurs when a and d are real numbers and either b 
or bi is also real. In these cases one real eigenvalue is Ai = — (|a — (ip + |6p)/2, 
and the corresponding eigenmatrix ai is, in case b is real, the selfadjoint matrix 

0"! = I ^ Q I . in the case b is imagmary it is cti = I ^ ^ 

• Quitting the boundary, starting from the lazy state p(0) = [ \ ) : The 



eigenvalue e{t) to the density matrix p(t), starting as e{t) = 0, can be analyzed 
in its behavior as an analytic function. Its Taylor expansion gives, as has been 
stated and proved in Lemma fTTl e{t) = and e{t) = 0. Continuing the procedure 
of solving the eigenvalue equation term by term according to the order of gives: 

^e(t)|,=o = i|a|>|>-rfr- 
Also this third order derivative vanishes if a = d, whereas, in this case 

4.3 Mixing and demixing 

The mixing property of a state can change, in the course of time, in any direction 
for general h. The case of dissipation can start with a pure state and will end up in 
a mixed state. In case of decay, when h has only one eigenvector to the eigenvalue 
zero, each mixed state will decay and approach the pure decay product. The 
path in S can also start from the boundary and go through the maximally mixed 
state u, the following evolution will then act as purifying. This means increase 
of entropy, followed by its decrease. Such cases have been found numerically in 
|NTP07] and in |D04j . If such an up-and-down process appears in parallel fashion 
in two parts of the system, the phase relations between them can also show this 
non-monotonic behavior. 

One may imagine a fast decay procedure with up-down evolution of entropy 
in a subspace, followed by a slower procedure, going over to another subspace. 
What is the end product of the first decay is then decaying again. It is of course 
possible that these happenings repeat several times, giving several ups and downs 
of entropy. From the mathematicians point of view there is the question of 
forming functions by superpositions of a finite number of exponential functions. 
In reality these things happen in nuclear processes. 



5 Summary and Conclusion 

We consider the quantum dynamical semigroups representing time evolutions of 
open systems with finite dimensional Hilbert space. In this paper we concentrate 
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on semigroups with a simple GKS-Lindblad generator. An inequality concerning 
the decay of eigenvalues is recognized as an important clue to the behavior of 
mixed states near the boundary of the set of states; the rank can not decrease in 
finite time. The complementary clue is a study on the increase of the rank of p{t), 
eigenvalues which start at zero getting positive. There are few exceptions to an 
increase in first order ~ t. For this analysis the characterization of structures is 
essential, representing the operator h which is used to build the generator Vh as 
an upper triangular matrix. In relation to the subspaces appearing as eigenspaces 
to h we define, in one to one correspondence, certain faces of S, the convex set 
of states, as "lazy" faces, since the rank of the p in such a face does not increase 
in first order ~ t. 

The three basic ways of non-invertible time evolution in physics, i.e. dis- 
sipation, decay and decoherence/dephasing are mirrored in the mathematical 
structures of h and of the generator V^, and the decomposition of h is the basis 
for a complete characterization of all the stationary states, including a concrete 
formula. We have moreover characterized the dephasing as not allowing circular 
paths of states, and we have characterized the stationary states as attractive. 
With a view on applications in "quantum engineering" we stated a possibility of 
inverting the task: find an evolution which leads to a given state. This task can 
be fulfilled for each state with one exception: The maximally mixed state. 

Hereupon we demonstrated the findings on two-dimensional systems. 

6 Appendix 

6.1 Complete positivity 

Dynamical maps of S into S have to preserve the trace. The dual maps are then 
"unital" i.e. they map 1 to 1. 

Stinespring's theorem states that cj) mapping B{Tt) to B{TC) is completely 
positive and unital iff 



The predual map, acting on a state represented by the density matrix p G BiH) 
is therefore 



From a physicists point of view, (HOj) or fj4T]) are the relevant equations for a 
quantum dynamical map, where B{T-C) is the algebra of observables for a system, 
and p2 G 8(7^2) is the state of the "reservoir". So we actually don't need the 
mathematicians definition of complete positivity. In the sequel we use only (HOll 
instead. 



3n2, 3U = iU^)-^onn^n2, 3p2eS(H2) 
VF: <I>(F) = Tr2((l®p2)-f/^-(i^®l)-f/). 



(40) 




(41) 
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18 PROPOSITION. The Kadison inequality. // $ is a unital completely 
positive map BiTi) BiH), then for each F G B{T-C) there holds the inequality 

<I>{F^F) > <I>{F^)<I>{F) (42) 

Proof. By doubling the "reservoir" -space 0.2 to 0.2 ® 'H2 one can consider the 
state p in B{T-C2) as the restriction of an entangled pure state, with a state vector 
in H2 ®'H2- Letting U act on H <^ H2 ® 'H2 one can reformulate (HO]) as 

= (^|f/"f • (F® 1 ® 1) • [/|^). (43) 

Now using 1^1') (^^l < 1 in $(^"1") • $(-F) gives the Kadison inequality. □ 

If one feels uneasy with the abstract algebraic notation, one may read (H3l) as 

((/)|$(F)|V') = {(j) ^\U^ ■ {F ® 1 ® 1) ■ U \ij ® (44) 

and use a set of basis vectors in H to write the product ^{F^) ■ $(F) in 
matrix notation, and insert |0„ (S> ^)(0n ® \E'| < 1 between ^{F^^) and $(-F). 



6.2 Block matrices 

Consider a splitting of H into subspaces, TC = TtiQ)Ti.2- Write vectors ip = X®' 
as ( ^ ) and write operators acting in Ti, as block matrices: 



A B 
C D 

R Q 

S 



(45) 
(46) 



The action of P/, : p ^ p=(^^^ 9^ (47) 

We need two special cases: 

The first one is the case of density matrices p with rank less than the dimension 
of the Hilbert space: 

If Q = and ^ = 0, then S = C RCK (48) 
The second special case with no restriction on p uses the Schur-Toeplitz triangu- 

lation, |L69j . of the generating matrix h: 

If C = 0, then 

R = Va{R)+BSB^ + AQB^ + BQ^A^ -llQB^A + A^BQ^ (49) 
Q = AQD^ + BSD^ - ^[A^AQ + QB^B + QD^D + RA^B + A^BS] (50) 
S = Vd{S)-1[B^BS + SB^B + B^AQ + Q^A^B]. (51) 
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For a general matrix 

( R Q\ 

which is in general not self adjoint, we have to replace in fl49l) and in flSip by 
M, and we also need 

M = DMA^ + DSB^ - I [MA^A + B^BM + D^DM + B^AR + SB^A] (52) 
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